The Einstein constraints: uniqueness and non-uniqueness in the conformal thin 

sandwich approach 
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We study the appearance of multiple solutions to certain decompositions of Einstein's constraint 
equations. Pfeiffer and York recently reported the existence of two branches of solutions for iden- 
tical background data in the extended conformal thin-sandwich decomposition. We show that the 
Hamiltonian constraint alone, when expressed in a certain way, admits two branches of solutions 
with properties very similar to those found by Pfeiffer and York. We construct these two branches 
analytically for a constant-density star in spherical symmetry, but argue that this behavior is more 
general. In the case of the Hamiltonian constraint this non-uniqueness is well known to be related to 
the sign of one particular term, and we argue that the extended conformal thin-sandwich equations 
contain a similar term that causes the breakdown of uniqueness. 
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I. INTRODUCTION 

With the help of a 3+1 decomposition Einstein's equa- 
tions can be split into a set of constraint equations and 
a set of evolution equations [1, 2]. The four constraint 
equations - one in the Hamiltonian constraint and three 
in the momentum constraint - constrain the induced spa- 
tial metric gij and the extrinsic curvature Kij on spatial 
hypersurfaces representing instants of constant coordi- 
nate time t. The constraint equations constrain only four 
of these initial variables; the remaining ones are freely 
specifiable and have to be chosen independently before 
the constraint equations can be solved. A decomposition 
of the initial data separates the freely specifiable vari- 
ables from the constrained ones. Given a particular de- 
composition, the construction of initial data then entails 
making well-motivated choices for the freely specifiable 
independent background data and then solving the con- 
straint equations for the constrained variables. 

The conformal thin-sandwich decomposition has 
emerged as a particularly popular decomposition among 
numerical relativists, especially for the construction of 
quasi-equilibrium data (see, e.g., the reviews [3, 4] and 
a brief discussion below). This variation of the original 
(non-conformal) thin-sandwich decomposition [5-7] was 
formally developed by York [8] (see also [9]), but more 
restricted versions had been introduced earlier [10, 11]. 
It has been used, for example, for the construction of bi- 
nary neutron stars (e.g. [12-15]), binary black holes (e.g. 
[16-20]) and black hole-neutron star binaries ([21, 22]). 

Given this wealth of experience with the conformal 
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thin-sandwich decomposition, it came as quite a sur- 
prise when Pfeiffer and York ([23], hereafter PY) re- 
cently discovered non-uniqueness in the solution of the 
conformal thin-sandwich system. Even for 'small' inde- 
pendent background data, which one would expect to 
generate gravitational initial data close to a flat slice of 
flat spacetime, the so-called "extended" set of confor- 
mal thin-sandwich data allowed for two branches of solu- 
tions. One of these two branches has comparatively weak 
gravitational fields and indeed approaches flat space in 
the limit of vanishing background data, while the second 
strong-field branch approaches a singular solution. 

In this paper we discuss some of the aspects of this 
non-uniqueness. We consider a spherically symmetric, 
constant density star and show that even for this very 
simple model we can find - for sufficiently small values 
of the density - two branches of solutions. These two 
branches share some of the characteristics of the solu- 
tions found by PY and illustrate their properties. For our 
spherically symmetric solution the non-uniqueness of so- 
lutions is caused by a particular term having the "wrong 
sign" (see, e.g. [2]), and we argue that the non-uniqueness 
found in the extended conformal thin-sandwich equations 
may be caused by a similar term. 

We note also that certain constrained evolution 
schemes [24, 25] solve a set of elliptic equations at ev- 
ery timestep which is very similar to the extended con- 
formal thin sandwich equations. These authors observed 
occasional failure of their elliptic solvers in the strong 
field regime, and it was argued [25, 26] that this fail- 
ure is caused by the "wrong sign" in the maximal slicing 
condition. 

This paper is organized as follows. In Section II we 
briefly review the conformal thin-sandwich decomposi- 
tion. Wc then consider the Hamiltonian constraint in 
spherical symmetry in Section III. First, in Sec. Ill A, 
we construct analytic solutions for constant density stars 
and show that these solutions consist of two branches 
with properties very similar to the solutions found by 
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PY. Subsequently, in Sec. IIIB, we prove that at least 
some of these properties persist for arbitrary spherically 
symmetric solutions. Finally, we provide a brief summary 
and discussion in Section IV. 



II. THE CONFORMAL THIN-SANDWICH 
DECOMPOSITIONS 

Conformal decompositions of the constraint equations 
start with a conformal transformation of the spatial met- 
ric, gij = ip^ijiji where ip is the conformal factor and g^ 
the conformally related metric. The Hamiltonian con- 
straint then becomes an equation for the conformal fac- 
tor 

V 2 V ~ \m ~ + \il>- 7 & Ay + 2^p = 0. (1) 

Here V and R are the covariant derivative and the trace 
of the Ricci tensor associated with g^ , and the extrinsic 
curvature is decomposed into its trace K and the confor- 
mally related trace-free part A % i , 



K ij = ip~ w A^ + \g ll K. 



(2) 



For completeness we have also included the matter source 
p = n a n b T a }>, where n a is the normal on the spatial hy- 
persurface and T a j, the stress-energy tensor, and where 
summation is carried out over four spacetimc indices. 

The matter term as written in Eq. (1) has the defect [2] 
that its positive sign combined with the positive exponent 
of ip prevent use of the maximum principle to prove local 
uniqueness of solutions. Therefore, it is not immediately 
clear that solutions to Eq. (1) are unique (indeed, we 
show in Sec. Ill, that often they are not unique). This 
defect can be cured [2] by introduction of a conformally 
scaled matter density p = ip 8 p; taking p > as freely 
specifiable data, the matter term in Eq. (1) becomes 
2TTip~ 3 p. Because of the sign-change in the exponent, 
this term is now well-behaved and the maximum princi- 
ple is applicable. We will use Eq. (1) with the unsealed 
p as a toy example in Sec. Ill below. Besides that, we 
are only interested in vacuum space-times and therefore 
do not include matter terms in the rest of this Section. 

The conformal metric g^, meanwhile, is freely speci- 
fiable. In the conformal thin-sandwich decompositions, 



the time derivative of the conformal metric, 



d, 



't9ij 



is also considered freely specifiable. Using the evolution 
equation for the spatial metric we can relate uu to , 



(3) 



where the conformal (or densitized) lapse ./V is related to 
the lapse N by N — ip 6 N. Inserting this expression into 
the momentum constraint yields 
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=-ip 6 V l K 



j ( — -u 

\2N 



(4) 



where (L/3) y = 2V^) - 2/3g lj V k (3 k is the conformal 
longitudinal operator. 

There are two versions of the conformal thin sandwich 
approach. In the standard conformal thin sandwich equa- 
tions, one specifies (gij,Uij; K, N) and suitable matter- 
terms, if applicable. Given these background variables, 
Eq. (1) and (4) (together with (3)) can be solved for the 
conformal factor ip and the shift (3 l , which completes the 
set of initial data. 

For maximal slices, K — 0, Eqs. (1) and (4) decou- 
ple, so that Eq. (4) can be considered first. For any 
given strictly positive N, this equation is a linear elliptic 
equation so that the existence of a unique solution (3' 1 is 
guaranteed. This is a key motivation for the entire struc- 
ture and is discussed in [8, 9]. Eq. (1) - with zero matter 
density p - becomes the standard Lichnerowicz equation 
for the conformal factor [27] and again has a unique solu- 
tion as long as the base metric is in the positive Yamabe 
class. 

In the extended system one regards dtK instead of TV 
as freely specifiable. The lapse can then be solved for 
from the trace of the evolution equation for the extrinsic 
curvature, which often is written as 



V 2 (7VV> 7 ) - (Nip 7 ) 



= -iP 5 (d t K-(3 k d k K). 



(5) 



The independent background data now are 
(gij,Uij- : K,d t K) (and suitable matter terms, if ap- 
plicable) and we solve five coupled elliptic equations (1), 
(4) and (5) for the conformal factor ip, the shift (3 l and 
the lapse N. This extended system has become very 
popular in numerical relativity because the ability to 



set the time dcrivati 



and dtK to zero provides a 



means of constructing quasi-cquilibrium data. 

However, PY demonstrated that the extended confor- 
mal thin sandwich equations behave very differently from 
the standard set, even for K = = d t K. Specifically, 
they found two branches of solutions for the same choices 
of free data. 

We wish to point out that Eq. (5) is written in a mis- 
leading way. As written, it appears that the maximum 
principle can be used for Eq. (5). However, A % i contains 
the lapse itself, cf. Eq. (3); displaying this dependence 
explicitly results in 



>T2 K ^ 



-iP 5 (d t K-/3 k d k K). 



(6) 



The first line of this equation has the structure 

v 2 (iv^ 7 )-/(^ 7 ) _1 

with non-negative coefficient /. The sign of / combined 
with the negative exponent of (Nip 7 ) in the second term 
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prevents application of the maximum principle, as did 
the unsealed density term in the Hamiltonian constraint 

(!)• 

We believe that this term might very well be responsi- 
ble for the complex behavior exhibited by the extended 
conformal thin sandwich equations. To support our 
claim, we analyze the Hamiltonian constraint (1) with 
an unsealed density in Section III below. We construct 
an analytic solution in spherical symmetry and explic- 
itly show the existence of two branches of solutions with 
properties very similar to those reported by PY. 



III. HAMILTONIAN CONSTRAINT WITH 
UNSCALED MATTER DENSITY 

As we discussed above, the Hamiltonian constraint 
Eq. (1) with unsealed matter density is not amenable 
to the maximum principle, and it turns out to be inter- 
esting investigate consequences of this fact. We consider 
the initial value problem at a moment of time-symmetry, 
Kij = 0, so that the momentum constraint is satis- 
fied identically. Assuming further conformal flatness and 
spherical symmetry, the Hamiltonian constraint Eq. (1) 
reduces to 



V 2 V> + 2ttpV 5 = 0, 
with ip > and with boundary conditions 



or 



1, 



r = 0, 
r — > oo 



(7a) 

(7b) 
(7c) 



where V 2 = d 2 /dr 2 + (2/r)d/dr represents the fiatspace 
Laplacian, and we assume a density profile p(r) > 0. 



Considering the function Cu ai we find that this function 
satisfies Eq. (7a) for any choice of a, given that C — 
(27rp /3) _1/4 . Indeed any solution ip to Eqs. (7) must be 
of this form in the interior of the star: The function Cu & 
with a = C 2 [ip(0)]~ 2 has the same value and derivative 
as ip at the origin, and as we show in the next section, 
this implies that ip = Cu & throughout the interior of the 
star. 

In the exterior, the only solutions of the flat-space 
Laplace equation with asymptotic value unity are the 
functions (3/r + 1, for some parameter (3. Consequently, 
any solution to Eq. (7) must be a member of the family 
of functions 



Cu a (r), r < R 

V>M = { , 1 ^ p 

— hi, r > R 



(11) 



with C given above, and a, (3 real parameters. The pa- 
rameters a and (3 are determined by continuity of ip and 
its first derivative at the surface of the star, 



| + 1 - Cu a (R), 
-jp=Cu' a {R), 



(12) 
(13) 



where a prime denotes d/dr. Eliminating (3, we find that 
a has to satisfy, 



PoR 2 = ^/», 



where 



/(«) = 



Or 



(1 + Q 



2N3 ' 



(14) 



(15) 



A. The constant density star 

First we will consider a constant density star of (con- 
formal) radius R and mass-density 



p(r) = 



Po, r < R, 
r>R. 



0. 



(8) 



We will take R to be fixed, and examine the solutions 
of this equation as we vary po- Thus, po plays the role 
of the "amplitude" of the perturbation away from trivial 
initial data. 

Solutions of Eq. (7) in the interior of the star can be 
found with the help of the so-called Sobolev functions 



u a (r) = 



{aR) 1 ' 2 



[r 2 + (aR) 2 



1/2 ' 



which satisfy 



V 2 u Q 



-3i£. 



(9) 



(10) 



Given a value for a we can find (3 from (12) or (13), which 
then completely specifies a solution to (7). 

The non-uniqueness of the solutions arises through the 
properties of the function f(a). We can see immediately 
that f{a) approaches zero for both a — > and a — > oo. 
For a sufficiently small value of poR 2 in (14) we may 
therefore pick either a small or a large value of a, which, 
as we will show below, corresponds to either a strong-field 
or a weak-field solution. 

Examining /(a) more carefully, we see that it takes 
its maximum at a c = \/5. Therefore, Eq. (14) has no 
solution if p is larger than the critical value 



Pc = 



2irR 2 



/ 2 (« c ) = 



3 5 5 
2irR 2 ¥ 



0.0320 

— W ' 



(16) 



At the critical density Eq. (14) has exactly one solution, 
a = ol c , while below the critical density there are two 
solutions; one with a < ce c and one with a > ce c . This 
behavior is in complete analogy to the behavior of the 
extended conformal thin-sandwich system examined in 
PY. 
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FIG. 1: ADM-energy and rest mass as function of poR 2 for 
the constant density star of Section III A. 



Having just derived all solutions to Eq. (7), we now 
now discuss their properties in more detail. It turns out 
to be convenient to parametrize these solutions by a. 
Each value of a corresponds to precisely one solution 
with po given by Eq. (14). Both limiting cases, a — > 
and a — > oo correspond to the limit of vanishing mass- 
density (see Eq. (14)). 

We begin by computing the ADM-energy, which can 
be found using Eq. (12), 



E = 2(3 



(17) 



For large a, the ADM-energy tends to zero and we re- 
cover flat space. In the limit a — > 0, however, the energy 
grows without bound, despite the fact that p — ► as 
a — ► 0. This establishes the a > a c branch as the weak- 
field branch, and a < a c as the strong-field branch. We 
show a graph of the energy as a function of density in 
Figure 1. 

Next we consider the rest mass M of the star, which 
is given by 

M= [ PoyfgdV = [ p i/j 6 4?rr 2 dr 

Jr<R JO 
3 

= — =■ la - a 5 + (1 + a 2 ) 3 arctan(a _1 )l R, (18) 
4a° 1 J 

where we have used Eq. (14) to eliminate po. This ex- 
pression has the limiting values 



Mk—R 

Of* 



for a — ► 0, 
for a, — > oo. 



(19) 
(20) 



The weak-field limit a — > oo corresponds to the limit in 
which the star has vanishing mass, whereas the strong- 
field limit a — > results in a star with unbounded mass, 
even though the density itself approaches zero. This be- 
havior is caused by the fact that the conformal factor, 
and hence the proper volume inside the stellar radius R 
diverges more rapidly than the rate at which the density 
po vanishes. 

We point out that for all a > we have E < M, so 
that the star has negative binding energy as expected. In 
the Newtonian limit a — > oo we recover the Newtonian 
binding energy, 



E 



M w — ^jR : 
5a 4 



3 M 2 
~5~R' 



(21) 



where we have used Eq. (20) in the second step. 

Finally, we locate the apparent horizons in this family 
of initial data sets. For a time-symmetric hypersurfacc, 
apparent horizons coincide with maximal surfaces, which 
in spherical symmetry and conformal flatness are given 
by the roots of 



9-0 tp 
dr 2r 



= 0. 



(22) 



For a > 1, no roots to this equation exist, so that the 
initial data surface does not contain an apparent hori- 
zon. For a < 1, two roots exist, one in the interior of 
the star at r = aR, and one in the exterior at r — R/a 2 . 
The latter one is the outermost extremal surface, which 
is the apparent horizon. Both extremal surfaces merge 
on the surface of the star for a = 1. The density, ADM- 
energy and rest mass at formation of the apparent hori- 
zon are p \ a =i = 3/(1287ri? 2 ) w 0.0075iT 2 , E = 2R, 
and M\ a=1 = 3nR/2 w 4.71i?. 

We now turn our attention to the critical point. 
Around its maximum f(a) behaves like a parabola, there- 
fore 



a - a c cx ±(p c - p ) 1/2 . 
Energy and mass at the critical point are 
2 



E r 



M r 



18 / 



(23) 



(24) 



125 



f9V5arctan(l/V5)-5ji?« 0.499i?, (25) 



respectively. Since dE/da ^ there, the energy also 
changes parabolically with p , 



E - E c cx a - a c cx ±(p c - p ) 



1/2 



(26) 



This parabolic behavior is apparent in Fig. 1. 

At the critical point the local uniqueness of solutions 
must break down, since the two branches meet there. For 
this to happen, the linearized operator must have a non- 
trivial solution at the critical point. The linearization of 
Eq. (7) reads 



V 2 5ip + 107rp i/> 4 5tp = 0, 



(27a) 



■5 



with boundary conditions 



0, 



dSip 
dr 
Sip^Q 



r = 
r — > oo. 



(27b) 
(27c) 



We will now construct all solutions of Eq. (27). While 
doing so, we consider p as given and fixed. If Sip = 
is the only solution, then the kernel of this equation 
is trivial, and solutions to the non-linear equation (7) 
are locally unique. As just argued, at the critical point 
this will not be the case, and there must be a non-zero 
solution of Eq. (27). As it turns out, we can construct 
this solution analytically. 

The key to solving Eq. (27) are again the Sobolev func- 
tions u a . Recall that ip — Cu a satisfies Eq. (7) in the 
interior of the star for any value of a. We can therefore 
take the derivative of Eq. (7) with respect to a and find 



2 du a 4 4 du a 

V z — h lOirpoC^ul— — = 0. 

da da 



(28) 



Choosing a to be a solution of (14), so that it is consistent 
with po, we can identify C A u\ — ip 4 , and Eq. (28) reduces 
to Eq. (27a). Consequently, any function Adu a /da, with 
a given by (14) and A an arbitrary constant, satisfies 
Eq. (27a). This forms a one-parameter family of func- 
tions, all of which automatically satisfy the differential 
equation Eqs. (27a) and the boundary condition (27b) in 
the interior. 

Solutions Sip in the exterior must satisfy the Laplace 
equation and the outer boundary condition (27c), i.e. 
they must take the form B/r for some constant B. Since 
this is a one-parameter family of solutions, we have found 
all solutions to Eqs. (27a) and (27c) in the exterior. 

To find a global solution Sip we now have to find con- 
stants A and B so that the interior solution Adu a /da 
matches the exterior solution B/r continuously in both 
the functions and their first derivatives at the stellar ra- 
dius r — R. As expected, non-trivial solutions with non- 
zero A and B exist only at the critical point a = a c . 
There, the solution Sip takes the form 



Sip c (r) oc < 



5R 2 - r 2 
(r 2 + 5i? 2 ) 3 / 2 : 
4 



I 6 3 / 2 r 



r <R 



> R 



(29) 



and it is easy to verify that it indeed satisfies Eq. (27) at 
the critical point. 



B. Results for general p > 

A fully worked out example like the constant density 
star presented above is very instructive. However, the ex- 
ample itself does not provide any indication whether its 
behavior is generic. In this Section we prove theorems 
valid for general p > with compact support, indicating 



that the behavior found for the constant density star is 
indeed generic for the Hamiltonian constraint with un- 
sealed matter density. We will first show that for suffi- 
ciently "large" matter-densities p, no solution exists. We 
will then consider the critical point and show that if a 
critical point exists, the solution must vary parabolically 
close to it, as did the constant density star, cf. Eq. (23). 
Finally, we will prove a result which was stated above to 
show that all solutions to the constant density star have 
been found: If two functions each satisfying Eqs. (7a) 
and (7b) have the same value at the origin, then they are 
identical. 

We start with some preliminaries. Rewriting the 
Laplacian in Eq. (7) we find 



(ip + rip')' = -2irrpip 5 < 0. 



(30) 



Therefore the combination ip + rip 1 is monotonically de- 
creasing and bounded from below by its asymptotic value 
for large r, 



ip + rip' > 1. 



(31) 



Furthermore, integrating Eq. (7) over a sphere of radius 
R we find 

AttR 2 ip'(R) = -2tt f pip 5 dV. (32) 
Jr 

Since p > and ip > we have ip' < 0, so that ip is 
a decreasing function of radius, which is bounded from 
below by its asymptotic value, ip > 1. 

We can now show that for sufficiently large p Eq. (7) 
does not admit strictly positive solutions ip. Solving 
Eq. (32) for ip' {R) and substituting into Eq. (31) we find 

2R[tp(R)-i\ > [ pif dV > ip{Rf~ k / ip k pdV, (33) 
Jr Jr 

for any k < 5, where the last inequality follows from 
ip' < 0. Rearranging terms we obtain 



(34) 



The right hand side of this inequality is bounded inde- 
pendently of the value of ip(R) by the biggest value of 
the function 



x - 1 

for x > 1. For k < 4, this function is bounded by 
so that any solution ip satisfies the integral bounds 



^J^ k pdV<C k 



(35) 



(36) 



(37) 
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for any k < 4. For a given k, the bound Ck is independent 
of R. For positive k, < k < 4, this inequality constrains 
how large solutions can be. For example, 



d > 



rIr 



ip A pdV > 



R 



[ pdV 

JR 



implies 



MR) < a 



1/4 



IIS/ 



dV 



-1/4 



(38) 



(39) 



which is a bound of how quickly the "upper" branch can 
diverge as p —* 0. 

For k = 0, the inequality (37) becomes independent of 
ip: If a solution ip exists for a certain p, then 



RJr 



0.163. 



(40) 



Equation (40) holds for any R for any strictly positive 
solution of Eq. (7), therefore if a density distribution p(r) 
satisfies 



m{R) 



f 4 4 
J pdV> 2—J; 



5 5 ' 



(41) 



even for one R, then no regular solution to the Hamilto- 
nian constraint Eq. (7) exists for this density. 

For the constant density star, m(r) is largest at the 
surface of the star, r = R, where 



m(R) 



l> {) I 4:irr 2 dr = ^-p R 3 - 
lR 3 



(42) 



Equation (41) then gives the necessary bound p < 
0.0389/i? 2 for the existence of solutions. Comparison 
with the exact critical density p c — 0.0320/i? 2 from 
Eq. (16) reveals that the upper bound of the theorem 
is only 20 per cent larger than the exact critical density 
(see also Fig. 1). 

Let us now examine the character of the critical point. 
We take a smooth sequence of non- negative densities, p 7 , 
such that p = when 7 = 0. We then look for a smooth 
sequence of solutions ip-y to Eq. (7) with the density p 
given by p~ n starting from ip = 1 at 7 = 0. The Implicit 
Function Theorem tells us that as long as the linearized 
equation, Eq. (27), 



v 2 Sip + 1071-^1/4 H = 0, 



with boundary conditions 



0, 



d&tp 
dr 



r = 



(43a) 

(43b) 
(43c) 



has no non-trivial solution, then the full nonlinear equa- 
tion 



V 2 V 7 + 27rp 7 ^ = 0, 



(44a) 



with boundary conditions 

7T = °> r = °' 
or 

tjjj — > 1, r — > 00, 



(44b) 
(44c) 



has a regular solution which changes smoothly as a func- 
tion of 7. 

The obvious question to ask is what happens if the 
sequence approaches the point where the first kernel of 
Eq. (43) appears 1 . Let us assume this happens at 70. 
The trick is to consider the limiting process rather than 
the limit point itself. We know that when 7 = 70 the 
equation 



v 2 e + ion Plo ^ o e 



0, 



(45) 



has a positive solution 9, going to zero at infinity. This 
is the ground state of a Schrodinger equation, because it 
is the first appearance of a kernel, hence it has no nodes, 
and thus it can be chosen to be everywhere positive. 

We now differentiate Eq. (44a) with respect to 7 (at 
any 7 < 70 to find 



2 #. 



^L + lOnp^ — 



4#7 = _ 2?r ^B 



dj 



(46) 



Multiplying Eq. (45) by # 7 /d7, Eq. (46) by 6 and sub- 
tracting the results we obtain 



v. ( ^vtf-ev^) = 

c?7 c?7 



dip. 



• dip. 



d~f d~f 



=10tt(p 7 V> 7 - P7oV4o)-^' 



2n^t 

d7 ' 



(47) 



Next we wish to integrate this equation over the whole 
space. Let us assume that p 1 has compact support, or, 
at least, falls off rapidly at infinity. Both 9 and dipj/dj 
fall off at infinity like 1/r and their first derivatives fall 
off like 1/r 2 . This means that the total divergence upon 
integration becomes a surface term, and the integrand 
falls off like 1/r 3 . Therefore the integral vanishes. We 
get 

5j( P ^-p^ )^9dV+ J d ^9dV = Q. (48) 

Consider this in the limit as 7 — > 70. The second term 
tends to the constant 



dp 1 
c?7 



(49) 



70 



Note that both 9 and ip la depend only on p~ (0 , and not 
its derivative. Therefore, changing dp 7 /d"f\ via 

Pi -> Pi + v{l - 70) (50) 



1 Clearly there are sequences p 7 for which this never happens, e.g. 
p 7 = 0. 
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for any function v(r) will change X as 

!->! + J vrfjdV. (51) 

Clearly, except for special instances, X will be non-zero. 
Let us now assume this generic case, 1^0. In the limit 
7 — ► 7o, the second term in Eq. (48) becomes the non- 
zero constant I, whereas the first term seems to go to 
zero. This cannot be and thus we are forced to conclude 
that the limiting process as 7 — > 70 must be somehow 
singular. 

The only thing that can possibly go bad is that 
dtp^/dj — > 00. Not only has it to blow up, it must do so 
over an extended region. This is the only way that the 
integral, in the limit, can go to a nonzero value. Let us 
assume 

^«( 7 o-7) P ( 52 ) 

for some negative power p < 0. This implies 

V> 7 -V>7o«(7o-7) P+1 - (53) 
Consider now the term 

p 7 V>* - P7oV4 = (P-y - P70) ^ 7 (54) 
+P70 (V> 7 + + V> 7 V> 70 + ^ 7o ) (^7 ~ V> 70 ) 

in the first integral in Eq. (48). The first term on the 
right hand side scales as 70 — 7 as the critical point is 
approached, whereas the second one scales as (7 — 7o) p+1 - 
Because p+1 < 1, the second term dominates, and the 
full integrand scales as 

{p,^ - p 70 <) ^ cx (70 - 7) p+1 (7o - 1? (55) 

close to the critical point. In the generic case, the integral 
has to approach the finite, non-zero value X in the limit 
7 — > 70, which can only happen if p = —1/2. Therefore, 
the solution must vary parabolically, 

^7 -^70 « (7 -7o) 1/2 • (56) 

This is exactly the behavior we have seen in the constant 
density star model and also with what Pfciffer and York 
[23] observed. 

This parabolic nature of solutions near the critical 
point can be demonstrated explicitly even in the non 
spherically symmetric case by using Lyapunov-Schmidt 
techniques [28]. We conjecture that a second branch of 
solutions exists beyond the critical point as a consequence 
of this parabolic behavior. We have seen this explicitly 
for the constant density stars, and we again refer to [28] 
for a more general treatment. 

So far we have shown that there are distributions p(r) 
for which no solutions of Eq. (7) exist, and based on the 
parabolic nature of the solutions at the critical point we 



have conjectured that there are distributions for which 
exactly two solutions exist. We do not know whether 
this is generic. 

Having moved past the first critical point, an open 
question is whether another critical point is reached. Im- 
mediately past the first critical point, p does not change, 
but the conformal factor increases. In the language of the 
Schrodinger equation this means that the potential deep- 
ens and the zero-energy ground state becomes a bound 
state with negative energy. As one moves away from the 
critical point along the upper branch one is moving 'back' 
toward smaller 7, and so we expect p 7 to decrease while ip 
continues to increase. One could have that pip 5 increases 
enough that the first excited state appears with zero en- 
ergy, or that pip 5 decreases again so that the ground state 
becomes a zero energy state again. In either case, the sys- 
tem reaches another critical point and the solution curve 
may turn again. Alternatively, p-ytp 5 may be such that 
neither of these two cases happens and the solution con- 
tinues on all the way to 7 = 0. This last alternative oc- 
curs for the constant density star, as we have shown by 
explicit calculation; however, we do not know whether 
this behavior is generic. 

Finally, we show that if we have two positive solutions 
tpi and ip2 to Eq. (7) whose maxima agree, then they are 
identical. To prove this we first note that the maxima of 
both must occur at r = 0. The maximum principle tells 
us that there cannot be a positive minimum, therefore 
there can only be one maximum, and therefore it must 
occur at the origin. Therefore at r = 0, the two func- 
tions tpi and ip2 agree, their first derivatives both vanish, 
and the second derivatives are equal (from Eq.(7)). By 
differentiating Eq.(7), one can show that all the deriva- 
tives of the two functions agree at r = 0. If the functions 
were analytic, we were done. However, there is no rea- 
son to expect that this be true. We need a more subtle 
argument. 

Track tpi and ip 2 as they move out from the origin. 
If they remain the same all the way to infinity, we are 
done. Instead, let us assume that at some point ipi > ip2- 
Therefore we must encounter a region in which d(ipi — 
ip2)/dr > and ipi — ip2 > and inside this region ipi > 
ip2- Take a point in this region, call it Rq. Consider the 
equations satisfied by ipi and ip2, subtract one from the 
other and integrate over the ball of radius Rq. We get 



f V 2 (^ 1 -^ 2 )rfT/ + 27r f ptyl - $)dV = 0. (57) 

J Rq J Rq 



The Laplacian becomes a boundary term, which is posi- 
tive, because the gradient of the difference is positive at 
Ro, while the bulk term is also non- negative. This cannot 
be, so therefore the initial assumption that the functions 
are different must be incorrect. 
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IV. SUMMARY AND DISCUSSION 

In this paper we investigate the reason for non- 
uniqueness in the extended conformal thin sandwich 
equations [23]. We argue that a term with the "wrong 
sign" in the elliptic equation determining the lapse, 
Eq. (6), is the cause for non-uniqueness. The sign of 
this particular term is such that the maximum principle 
cannot be applied to prove local uniqueness of solutions. 
We support our claim by examining a simpler equation 
having a term with an analogous "wrong sign", namely 
the the Hamiltonian constraint with unsealed [2] matter 
source p (cf. Eq. (1)). Specializing to constant den- 
sity stars we construct analytical solutions. We find two 
branches of solutions - a weak-field and a strong-field 
branch - with properties that are remarkably similar to 
those found by PY. 

We comment briefly that solutions to the original con- 
formal thin-sandwich decomposition, consisting of the 
Hamiltonian constraint (1) and the momentum con- 
straint (4) only, are unique [8, 23, 29-32]. PY found 
multiple solutions only for the extended conformal thin- 
sandwich decomposition, which includes the lapse equa- 
tion (5) in addition to the two constraints. This, too, 
suggests that the non-uniqueness is caused by the lapse 
equation, in accordance with our findings. 

Our findings are certainly relevant for numerical work: 
If one wants to solve the extended conformal thin sand- 
wich equations, then apparently, the possibility of find- 
ing two solutions is unavoidable. Whether this will pose 
a problem for numerical work is less clear. Sufficiently 
far away from the critical point, the solutions along the 
upper and lower branch are significantly different, and it 
should be obvious which solution is desired (generally 
the "lower" one, which reduces to flat space for triv- 
ial free data). Many different researchers have solved 
the extended conformal thin sandwich equations without 



problems [12-22] and have obtained a solution with sat- 
isfactory properties. However, past success is no guar- 
antee for future success, and for choices of free data 
which may be interesting in the future, non-uniqueness 
issues could very well arise, especially if one is interested 
in solutions which happen to be "close" to the critical 
point. Indeed, in constrained evolutions schemes, which 
solve elliptic equations similar to the extended confor- 
mal thin sandwich equations, it was reported that the 
elliptic solver failed to converge in near-critical collapse 
of Brill waves [24, 25]. It was further argued that this 
failure related to the "wrong sign" in a term of the max- 
imum slicing condition [25, 26]. How precisely a numer- 
ical code behaves in such cases depends very sensitively 
on its implementation. Some algorithms may not con- 
verge at all, like the multi-grid schemes in [24, 25], while 
other algorithms may converge to one of the two solu- 
tions (e.g. the Newton-Raphson method used in PY; cf. 
[33]). In the former case it may be difficult to ascertain 
whether failure of the numerical method is indeed due 
to non-uniqueness properties of the underlying analytic 
problem (rather than a bug), whereas in the latter case 
one is faced with the question of which of the two solu- 
tions one wants, and how to ensure convergence toward 
the desired solution. 
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